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MODELING  RELATIVE  POSITION,  RELATIVE  VELOCITY, 
AND  RANGE  RATE  FOR  FORMATION  FLYING 

Craig  A.  McLaughlin,  Chris  Sabol,'  Aaron  Swank, + 

Richard  D.  Burns  and  K.  Kim  Luu^ 

The  relative  position,  relative  velocity,  and  range  rate  evolution  is  examined  for 
various  formations  of  satellites.  A  simple  analytical  model  including  Earth  oblate¬ 
ness  effects  for  the  equations  of  relative  motion  is  presented.  This  model  provides 
physical  insight  into  the  Earth  oblateness  effects  that  are  neglected  by  using  Hill’s 
equations.  The  accuracy  of  the  relative  position,  relative  velocity,  and  range  rate 
predictions  for  the  analytical  model  are  compared  to  realistic  force  modeling  ob¬ 
tained  using  the  Draper  Semianalytical  Satellite  Theory  for  formations  of  varying 
size,  inclination,  and  altitude. 


INTRODUCTION 

Several  misginns  flying  satellites  in  formation  are  under  development  or  have  been  recently 
launched.  The  relative  velocity  and  range  rate  have  lately  become  a  larger  concern  for  formation 
flying  missions.  How  and  Tillerson7  show  that  expected  relative  velocity  measurement  errors  of 
2-3  mm/s  have  a  significant  effect  on  the  relative  estimation  and  control  problem.  This  causes  an 
increase  in  the  amount  of  fuel  required  to  maintain  a  formation  of  satellites.  One  solution  to  this 
problem  is  to  provide  improved  relative  navigation  measurements  through  an  intersatellite  range 
and  range  rate  sensor.  Such  a  sensor  may  be  subject  to  constraints  on  the  allowable  range  and  range 
rate.  Most  analyses  have  concentrated  on  the  relative  position  accuracy.  This  paper  presents  an 
analysis  of  relative  position,  relative  velocity,  and  range  rate.  In  addition,  a  simple  analytical  model 
is  developed  that  provides  physical  insight  into  Earth  oblateness  effects  that  are  ignored  when  using 
Hill’s  equations.  This  model  was  developed  to  address  the  need  for  a  way  to  describe  the  relative 
position,  velocity,  and  range  rate  of  the  satellites  in  a  formation  that  attained  a  high  level  of  accuracy, 
but  was  easy  to  use  and  understand.  In  addition,  the  model  provides  physical  insight  into  the  effects 
of  J2  on  the  relative  equations  of  motion  using  classical  descriptions  of  orbital  mechanics. 

Sabol  et  al.12  presented  a  numerical  analysis  of  satellite  formation  flying  that  showed  the 
absolute  and  relative  effects  of  Ji  are  sigificant  for  formations  of  satellites  with  out-of-plane  motion. 
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In  addition,  they  presented  a  simple  analytical  description  of  the  reason  for  the  observed  motion  in 
terms  of  J 2-  However,  that  analysis  only  looked  at  polar  orbits.  Alffiend  et  al.1  developed  a  state 
transition  matrix  to  analyze  the  errors  in  using  Hill’s  equations  for  differences  in  orbital  elements. 
Gim  and  Alfriend5  develop  this  analysis  further  to  produce  a  state  transition  matrix  of  the  relative 
motion  including  eccentric  reference  orbits  and  the  effects  of  Ji-  Inalhan  and  How8  present  an 
analysis  using  a  closed  form  solution  for  relative  motion  in  eccentric  orbits  as  a  function  of  true 
anomaly.  They  rely  on  a  development  originally  attributed  to  Lawden.9 

This  paper  presents  an  analysis  of  the  relative  position,  relative  velocity,  and  range  rate  of 
satellites  in  low  eccentricity  orbits.  If  x,  y,  and  z  represent  the  distance  between  satellites  in  the 
radial,  along-track,  and  cross-track  directions  and  x,  y,  and  i  represent  the  components  of  relative 
velocity  in  those  directions  then  the  range  p,  relative  velocity  v,  and  the  range  rate  p  are 


(1) 

p2  =  x2  +  y2  +  z2 

(2) 

v 2  =  x2  +  y2  +  i2 

(3) 

p  =  (xx  +  yy  +  zz)fp 

These  equations  hold  regardless  of  the  method  used  to  calculate  the  relative  positions  and  velocities. 

In  general  the  differential  equations  describing  the  relative  motion  of  satellites  where  the  refer¬ 
ence  satellite  is  in  a  near  circular  orbit  and  the  relative  separation  is  small  compared  to  the  radius 
of  the  orbit  are 

(4)  x  -  2wxvy  -  3u)2„x  =  fx 

(5)  y  +  2wxvx  =  /„ 

(6)  z  +  wjz  =  }z 

where  wxy  is  the  frequency  of  the  in-plane  motion,  wz  is  the  frequency  of  the  out-of-plane  motion,  and 
fx,  /„,  and  fz  represent  the  relative  perturbing  forces.  For  a  spherical  Earth  as  in  Hill’s  equations, 
wxv  —  wz  —  n,  where  n  is  the  mean  motion. 

The  Draper  Semianalytical  Satellite  Theory  was  used  to  model  the  truth  of  the  mean  element 
evolution  of  the  satellite  formations  using  realistic  dynamics.  An  analytical  model  was  developed 
that  accurately  captures  the  effects  of  J2  on  the  formation  evolution.  Results  of  the  model  were 
compared  to  DSST  for  formations  with  a  variety  of  inclinations.  In  addition,  the  effects  of  altitude 
and  formation  size  are  examined. 


Hill’s  Equations 


Hill’s  equations6  are  commonly  used  for  formation  design  and  analysis.  Sabol  et  al.,12  Alfriend 
et  al.,1  and  Inalhan  and  How8  point  out  some  of  the  limitations  of  using  Hill’s  equations.  These 
limitations  arise  from  assuming  the  reference  orbit  is  circular  and  the  Earth  is  a  perfect  sphere. 


Vallado16  provides  a  detailed  derivation  of  Hill’s  equations  and  the  solution  for  unperturbed 
motion  for  both  relative  position  and  velocity.  The  solutions  for  unperturbed  motion  are 


(7) 

(8) 
0) 
(10) 
(11) 
(12) 


x  =  (io/n)  sinnt  -  [3xo  +  (2y0/n)]  cos  nt  +  4x0  +  2 y0/n 

y  =  (2io/n)  cos  nt  +  [6x0  +  (4y0/n)]  sin  nt  -  (6nxo  +  3 y0)t  -  2x0/n  +  yo 

z  =  (i0/n)  sin  nt  +  z0  cos  nt 

x  =  x0  cos  nt  +  (3nx0  +  2y0)  sin  nt 

y  =  -2x0  sinnt  +  (6nx0  +  4y0)  cos  nt  -  (6 nx0  +  3yo) 

i  =  i0  cos  nt  -  nzo  sin  nt 
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Table  1:  Initial  conditions  for  circular  and  projected  circular  formations 


Formation 

xo 

X 0 

yo 

ifo 

zo 

zo 

Circular 

Projected  Circular 

(po/2)cosfl 
(Po/ 2)cos0 

-(pon/2)  sin0 
-(pf,n/2)sin0 

2  ±0/n 

2  ±o/n 

—2nxo 

—2nxo 

±\/3a:o 

±2xq 

±\/3£q 

±2±o 

where  the  subscript  0  values  are  initial  conditions.  Vallado16  and  Sabol  et  al.12  provide  more  detailed 
analyses  of  the  motion  that  results  from  Hill’s  equations.  Of  particular  interest  for  this  study  are 
the  circular  and  projected  circular  formations  presented  in  Sabol  et  al.12  The  circular  formation  has 
relative  motion  in  a  circle  about  the  reference  trajectory.  The  projected  circular  formation  projects 
a  circle  on  the  yz-plane  and  thus  onto  the  Earth  below.  The  initial  conditions  for  these  formations 
cure  given  in  Table  1.  po  is  the  radius  of  the  circle,  p'0  is  the  radius  of  the  projected  circle,  and  6  is  the 
phase  angle  in  the  circle  or  projected  circle  measured  counterclockwise  from  the  z  axis.  The  yo  and 
yo  conditions  set  the  along-track  offset  and  drift  to  zero.  The  signs  on  the  zo  and  iq  conditions  must 
be  the  same.  There  are  two  planes  in  which  both  the  circular  and  projected  circular  formation  are 
possible.  Both  intersect  the  cross-track/along-track  plane  along  the  along-track  axis.  The  circular 
formation  is  inclined  to  that  plane  at  ±30°.  The  projected  circular  is  inclined  at  ±26.565°  to  that 
plane. 


RELATIVE  MOTION  INCLUDING  J2 


The  primary  perturbation  affecting  formation  flying  satellites  of  identical  area-to-mass  ratio  is 
Earth  oblateness.  This  is  represented  by  the  factor  J?.  The  first  order  secular  effects  of  Ji  on  the 
right  ascension  of  the  ascending  node  fl,  the  argument  of  perigee  w,  and  the  mean  anomaly  M  are 


given  by 

(13) 

si 

(14) 

u> 

(15) 

M 

cos* 


3nR%J2 

32SA(4-5»inJi) 
v  v  ’ 


where  Re  is  the  equatorial  radius  of  Earth,  e  is  the  eccentricity,  p  is  the  semilatus  rectum,  and  i 
is  the  inclination.  As  reported  in  Sabol  et  al.12  the  J2  perturbations  have  two  primary  effects  on 
satellites  flying  in  formation.  The  first  is  a  long  periodic  effect  because  the  in-plane  (xy)  motion 
has  a  period  that  depends  on  perigee  passing,  but  the  out-of-plane  (z)  motion  has  a  period  that 
depends  on  nodal  crossing.  The  second  effect  applies  to  satellites  with  inclination  differences.  In 
this  case  the  above  equations  show  that  there  are  differential  perturbations  which  cause  secular 
growth  in  the  separation  of  the  satellites.  The  remainder  of  this  section  develops  successively  more 
complex  equations  to  describe  the  relative  motion  of  satellites  affected  by  Ji-  The  limitations  that 
the  satellites  are  in  near  circular  orbits  and  that  the  relative  separation  between  the  satellites  is 
small  relative  to  the  orbit  radii  will  still  apply  to  these  equations. 


Satellites  With  Matching  Inclinations 

Differences  in  either  right  ascension  of  the  ascending  node  fl  or  inclination  i  provide  out-of- 
plane  (z)  motion  described  in  Hill’s  equations.  In  both  cases  the  following  effects  will  be  present. 
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According  to  Hill’s  equations  all  three  components  have  motion  with  a  frequency  governed  by  the 
mean  motion  n.  However,  the  frequency  of  the  x  and  y  components  are  actually  determined  by  the 
satellites  passing  perigee.  For  orbits  perturbed  by  J2  this  is  governed  by  (n  +  M).  The  frequency 
of  the  z  component  is  actually  determined  by  nodal  crossing  time.  This  frequency  is  governed  by 
(n  +  M  +w)  for  orbits  perturbed  by  J2.  The  resulting  motion  is  described  by 


(16)  z 

(17)  y 

(18)  * 

(19)  x 

(20)  y 

(21)  i 


(±o/t uxy)  sin  wxyt  -  (3x0  +  2y0/wxv)  cos  wx1lt  +  4x0  +  2y0/wxy 

(6x0  +  4y0/wxy)  sin  wxyt  +  (2x0/wxy)  cost uxyt  -  (6wxyx0  +  3 y0)t  +  yo~  2x0/wxy 

Zocoswzt+  (zo/wx)smwzt 

io  costnZJ ,f  +  (Zwxyx0  +  2  y0)  sin  wxyt 

(6wxyx0  +  4y0)  cos  wxyt  -  2i0  sin  wxyt  -  (6wxyx0  +  3yo) 

-wzZq  sin  wzt  +  io  cos  wzt 


where 

(22)  wxy  =  n  +  M 

(23)  wz  —  n  +  M  +  Cj 

These  are  the  same  as  the  solutions  to  Hill’s  equations  except  that  wxy  replaces  n  in  the  x  and  y 
terms,  and  wz  replaces  n  in  the  z  term.  The  relative  velocities  are  simply  the  time  derivatives  of 
the  relative  positions.  The  initial  conditions  in  Hill’s  equations  were  used  to  determine  the  initial 
mean  orbital  elements  in  DSST.  However,  the  initial  conditions  in  the  model  are  calculated  using 
wxy  and  wz  instead  of  n  as  in  Table  1.  The  relative  separation  in  each  component  will  not  grow  with 
time,  but  the  phasing  of  the  motion  will  cause  the  overall  separation  to  have  a  long  period  effect. 
These  equations  predict  that  there  will  be  no  separation  growth  at  critical  incination.  For  satellites 
that  have  z  motion  caused  only  by  differences  in  ft  these  equations  accurately  describe  the  relative 
motion  in  all  three  components. 


Satellites  With  Inclination  Differences 


Differential  effects  on  ft,  ui,  and  M  occur  when  the  satellites  have  different  inclinations.  In  this 
case  J2  affects  both  the  y  and  z  motion.  In  the  z  equation  an  additional  term  appears  containing 
Silt.  This  term,  which  represents  the  effect  of  the  differential  precession  of  ft,  will  slowly  grow  with 
time.  However,  the  y  equation  has  an  additional  term  that  contains  all  the  differential  rates  (<Sft, 
5u,  and  SM)  multiplied  by  time.  The  <5ft  part  of  the  term  represents  the  differential  precession  of  ft. 
The  rest  of  the  term  represents  the  along-track  effects  of  differential  rates  of  oj  and  M.  This  term 
may  grow  quickly  with  time  depending  on  the  inclination.  The  final  set  of  equations  is  now 


(24)  x 

(25)  y 

(26)  2 

(27)  x 

(28)  y 

(29)  i 


(x0/wxy)smwxyt  -  (3x0  +  2y0/wxy) cos wxyt  +  4x0  +  2y0/wxy 

(6x0  +  4y0/wxy)sinwxyt  +  {2x0/wxy)coswxyt  -  {6wxyx0  +  3yo)t  +  yo  -  2x0/wxy 

+a(Sil  cos  i  +  Sih  +  8M)t 

zocoswzt  +  (io/wz)smwzt  —  aSilt  sin  i  cos  wzt 

i0  cos  wxyt  +  (3wxyx0  +  2y0)  sin  wxyt 

(6wxyx0  +  4y0)  cos  wxyt  -  2io  sin  wxyt  -  (6wxyx0  +  3 y0)  +  a(<5ft  cos  i  +  Su  +  SM) 
— w2zosiawzt  +  iocoswzt  —  aStl  sini  [cosu)z<  —  twz  sinu>2f] 


Because  of  the  perturbation  terms  the  initial  conditions  on  the  velocity  will  have  a  bias  compared 
to  DSST  in  the  y  and  z  terms  for  satellites  that  have  inclination  differences.  This  will  be  apparent 
in  the  results  presented  later. 
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Figure  1:  Perturbations  caused  by  J2  as  a  function  of  inclination  for  an  800  km  altitude,  1  km  radius 
circular  formation 


If  e2  terms  are  ignored  and  k  =  3nR2BJ2 /2p2  the  differential  rates  caused  by  J2  of  fi,  w,  and 
M  are 

(30)  Stl  =  k  sin  iSi 

(31)  6<jj  =  —5ksmicosiSi 

(32)  5M  =  — 3fc  sin  t  cos  j<5i 

The  secular  term  caused  by  J2  in  the  y  equation  then  becomes 

(33)  a(SQcosi  +Su>  +SM)t  =  -7 atk  sin  i  cos  iSi 

This  equation  shows  that  the  secular  term  in  y  is  zero  for  inclinations  of  0°,  90°,  and  180°.  The 
maximum  growth  will  occur  at  inclinations  of  45°  and  135°.  The  new  term  in  the  z  equation  is 
periodic  with  a  secularly  increasing  amplitude.  Substituting  in  this  term  yields 

(34)  —  a6(lt  sin  i  cos  wzt  =  —at  cos  wztk  sin2  iSi 

This  equation  shows  the  term  will  be  maximum  for  a  polar  orbit  and  zero  for  an  equatorial  orbit. 
The  differential  J2  terms  used  in  the  results  section  were  found  by  differencing  the  actual  rates  from 
Eqs.  13  through  15  for  each  orbit  instead  of  using  Eqs.  30  through  32.  Fig.  1  shows  the  relative 
velocity  magnitude  added  by  the  J2  perturbation  to  the  y,  z,  and  total  relative  motion  as  a  function 
of  inclination.  The  figure  shows  that  except  for  high  inclinations  the  y  perturbation  is  the  dominant 
perturbation.  Also,  the  maximum  perturbation  occurs  at  an  inclination  of  just  over  45°.  However, 
the  z  perturbation  is  out-of-plane  and,  therefore,  will  probably  be  more  expensive  to  correct  than 
the  y  perturbation. 
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PROPAGATIONS 


Test  cases  for  each  of  the  four  formation  designs  were  propagated  for  four  days  in  the  presence 
of  realistic  dynamics.  The  DSST  Averaged  Orbit  Generator  (AOG)  was  used  to  propagate  mean 
elements  for  each  of  the  test  cases.  DSST  was  chosen  for  its  rigorous  modeling  of  the  averaged 
equations  of  motion.  This  allows  for  fast,  efficient,  and  accurate  orbit  propagation  over  long  periods 
using  step  sizes  on  the  order  of  one  day.  McClain,10  Danielson,2  Feiger,3  Sabol  et  al.,13  and  Fonte4 
provide  more  information  on  the  formulation  and  accuracy  of  the  DSST.  The  DSST  runs  were  exe¬ 
cuted  with  the  Draper  Research  and  Development  version  of  the  Goddard  Trajectory  Determination 
System  (DGTDS).4  The  mean  elements  were  propagated  because  the  mean  evolution  of  the  motion 
is  what  is  important  for  comparing  the  satellites  in  a  formation.  Trying  to  control  the  periodic 
changes  in  the  relative  motion  caused  by  periodic  changes  in  the  osculating  orbited  elements  would 
not  be  advisable  unless  the  satellites  are  extremely  close  and  there  is  a  risk  of  collision.  However, 
Proulx  et  al.11  provide  an  examination  of  the  short-periodic  effects  on  formation  flying  of  satellites 
with  different  area-to-mass  ratios. 

The  perturbations  modeled  in  the  propagations  include  a  geopotential  to  degree  and  order 
21,  atmospheric  drag  effects,  luni-solar  third-body  point-mass  effects,  and  solar  radiation  pressure. 
Table  2  lists  the  specifics  of  the  ephemeris  generations. 


Table  2:  DSST  propagation  description 


Characteristic 

Setting 

Start  time 

Sept.  15,  1998,  0  hr,  0  min,  0  sec 

End  time 

Sept.  19,  1998,  0  hr,  0  min,  0  sec 

Integrator 

12th-order  summed  Cowell/ Adams  predict  partial  correct 

Step  size 

12  hrs 

Geopotential 

21x21  JGM2 

Atmospheric  drag 

Jacchia-Roberts,  Cd=2.0 
“Hot”  Schatten  Solar  Flux,  kp  prediction 

Third-body 

Solar/lunar  point  masses  based  on  JPL  DE  ephemerides 

Solar  radiation  pressure 

Cylindrical  shadow  model,  (7r=1.2 

The  Jacchia-Roberts  atmospheric  density  model  is  dependent  upon  daily  F10.7  cm  solar  flux  and 
3-hourly  values  of  the  kp  geomagnetic  indices.  Schatten14, 15  developed  a  method  to  predict  smoothed 
estimates  of  these  parameters  over  the  solar  cycle.  The  solar  flux  and  geomagnetic  indices  used  in  this 
analysis  are  the  “hot”  or  maximum  value  predictions  from  1994.  The  “hot”  atmosphere  prediction 
provides  a  conservative  estimate  of  drag  effects.  Identical  area-to-mass  ratios  of  0.01  m2/kg  were 
used  for  each  spacecraft  in  drag  and  solar  radiation  pressure  calculations. 

For  each  of  the  formation  flying  designs,  three  orbits  were  propagated:  one  representing  the 
reference  satellite,  and  the  other  two  representing  other  satellites  in  the  formation.  The  trajectories 
were  then  compared  over  four  days  to  examine  the  relative  motion  of  the  formations  in  the  presence 
of  perturbations.  The  output  of  the  relative  velocities  from  DGTDS  (xgtds>  Vgtds ,  zqtds )  was 
simply  the  difference  between  the  relative  velocity  components  of  the  satellites  in  a  intertial  frame. 
Consequently,  the  relative  velocities  had  to  be  adjusted  into  Hill’s  frame  by  subtracting  out  the  <3xr 
term.  The  rotation  rate  is  the  mean  motion  n  about  the  z  axis.  Therefore,  the  relative  velocities 
were  adjusted  by 

(35)  Amu  =  Adgtds  +  ny 


2170 


(36) 

(37) 


VHill 

zhui 


Vdgtds  -  nx 
zdgtds 


The  elements  resulting  from  the  formation  designs  and  input  into  the  DGTDS  DSST  propagator 
are  assumed  to  be  averaged  or  mean  elements.  If  the  Keplerian  elements  are  assumed  to  be  osculating 
in  the  presence  of  perturbations,  the  mean  semimajor  axes  will  not  match,  and  the  formation  design 
specifications  will  not  hold. 

The  basic  formations  studied  were  at  800  km  altitude  with  a  circle  or  projected  circle  with 
a  radius  of  1  km.  The  satellites  were  assumed  to  be  initially  at  the  ascending  node.  In  all  cases 
satellite  1  is  at  the  origin  of  the  Hill’s  frame,  satellite  2  is  at  9  =  270°,  and  satellite  3  is  at  9  =  90°. 
This  means  that  satellite  2  has  out-of-plane  motion  primarily  caused  by  inclination  differences  and 
satellite  3  has  out-of-plane  motion  primarily  caused  by  differences  in  right  ascension.  The  basic 
formations  were  examined  at  reference  inclinations  of  1°,  28.5°,  45°,  63.435°,  and  90°.  Each  satellite 
had  semimajor  axis  a  of  7178.1363  km.  For  all  formations  satellite  1  always  had  eccentricity  of 
lxlO-8,  reference  inclination,  and  zero  right  ascension  fl,  argument  of  perigee  w,  and  mean  anomaly 
M.  For  all  of  the  basic  formations  satellite  2  always  had  eccentricity  of  6.965597xl0~5,  and  satellite 
3  always  had  eccentricity  of  6.965463xl0-5  for  the  circular  formations  and  6.965705xl0~s  for  the 
projected  circular  formations.  The  other  initial  mean  orbital  elements  for  satellites  2  and  3  are  given 
in  Table  3. 

In  addition  to  the  basic  800  km  altitude,  1  km  circular  and  projected  circular  formations,  three 
additional  cases  were  examined.  The  first  was  a  1  km  projected  circular  formation  at  600  km  altitude 
(a  =  6978.1363)  and  28.5°  inclination.  The  other  two  were  projected  circular  formations  at  800  km 
altitude  and  28.5°  inclination.  One  projected  circle  had  a  radius  of  10  km  and  the  other  had  a  radius 
of  100  m.  All  satellites  in  a  given  formation  had  matching  mean  semimajor  axes.  The  initial  orbital 
elements  of  satellites  2  and  3  for  the  other  formations  are  given  in  Table  4. 


RESULTS 


The  results  presented  are  for  circular  formations  with  radii  of  1  km  at  800  km  altitude  unless 
specifically  stated  otherwise.  Satellites  1  and  3  have  out-of-plane  relative  motion  primarily  because  of 
differences  in  right  ascension.  Satellites  1  and  2  have  out-of-plane  relative  motion  primarily  because 
of  differences  in  inclination. 


Satellites  1  and  3 

Since  satellites  1  and  3  have  little  or  no  inclination  differences  compared  to  satellites  1  and  2, 
it  is  illustrative  to  examine  these  first  starting  with  a  polar  reference  orbit.  Figs.  2  through  5  show 
the  results  for  the  relative  separations.  Each  figure  shows  the  evolution  obtained  from  DSST,  the  J 2 
model,  and  the  error  of  the  J2  model  compared  to  DSST.  Although  the  total  separation  grows,  the 
separation  in  each  component  has  only  periodic  changes  and  no  secular  growth.  The  errors  of  the  J-i 
model  in  the  radial  component  are  less  than  a  meter,  and  the  errors  in  the  cross-track  component 
are  a  few  meters.  By  far  the  largest  error  is  in  the  along-track  component.  The  errors  shown  here 
are  typical  of  all  the  cases  rim. 

Figs.  6  through  9  show  the  results  for  the  relative  velocities.  Again,  each  figure  shows  the 
evolution  obtained  from  DSST,  the  J2  model,  and  the  error  of  the  J2  model  compared  to  DSST. 
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Table  3:  Initial  mean  orbital  elements  for  satellites  2  and  3  (Basic  formations) 


Formation 

*,  deg 

n,  deg 

u,  deg 

M,  deg 

Circ. 

t  =  1° 

Sat  2 

1.00691260 

5.480xl0_® 

270.02478675 

89.97515845 

Sat  3 

1.00002389 

0.39610432 

359.60395601 

0 

i  =  28.5° 

Sat  2 

28.50691260 

2.09xl0-« 

270.02483977 

89.97515845 

Sat  3 

28.50000077 

0.01448802 

359.98726768 

0 

i  =  45° 

Sat  2 

45.00691260 

1.48x10-® 

270.02484058 

89.97515845 

Sat  3 

45.00000042 

9.77658xl0-3 

359.99308692 

1.21x10-® 

t  =  63.435° 

Sat  2 

63.44191260 

1.21x10-® 

270.02484106 

89.97515845 

Sat  3 

63.43500021 

7.72906xl0-3 

359.99654347 

0 

O 

O 

02 

II 

*«* 

Sat  2 

90.00691260 

1.21x10-® 

270.02484155 

89.97515845 

Sat  3 

90 

6.91308xl0“3 

0 

1.21x10-® 

Proj.  Circ. 

*  =  1° 

Sat  2 

1.00798199 

6.321x10-® 

270.02478675 

89.97116746 

Sat  3 

1.00003186 

0.45737944 

359.54269022 

1.21x10-® 

t  =  28.5° 

Sat  2 

28.50798199 

2.41x10-® 

270.02883049 

89.97116746 

Sat  3 
i  =  45° 

28.50000102 

0.01672932 

359.98529799 

0 

Sat  2 

45.00798199 

1.71x10-® 

270.02883143 

89.97116746 

Sat  3 

45.00000056 

0.01128902 

359.99201746 

359.99999879 

i  —  63.435° 

Sat  2 

63.44298199 

1.21x10-® 

270.02883198 

89.97116746 

Sat  3 

63.43500028 

8.92475-3 

359.99600874 

0 

i  =  90° 

Sat  2 

90.00798199 

1.21x10-® 

270.02883254 

89.97116746 

Sat  3 

90 

7.98254xl0-3 

0 

0 

Table  4:  Initial  mean  orbital  elements  for  satellites  2  and  3  (Other  formations) 


Formation 

e 

i,  deg 

fl,  deg 

u>,  deg 

M,  deg 

Proj.  Circ. 
i  =  28.5° 

Alt. =600  km 
Sat  2 

7.165237x10-® 

28.50821076 

2.56x10-® 

270.02983170 

89.97016613 

Sat  3 

7.165374x10-® 

28.50000108 

0.01720883 

359.98487658 

0 

p'0  =  10  km 
Sat  2 

6.9656638xl0-4 

28.57981973 

2.3245xl0-4 

270.31876420 

89.68103165 

Sat  3 

6.9704276xl0-4 

28.50010254 

0.16739767 

359.85288811 

0 

p'0  =  100  m 
Sat  2 

6.96560x10-® 

28.50079820 

0 

269.97223875 

90.02776123 

Sat  3 

6.96188x10-® 

28.50000001 

1.67283xl0-3 

359.99852989 

0 
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Figure  2:  Comparison  of  the  radial  separation  (x)  of  satellite  1  and  3  from  DSST  and  J i  model  for 
a  800  km,  polar  reference  orbit.  The  satellites  are  in  a  1  km  circular  formation. 


Figure  3:  Comparison  of  the  along-track  separation  (y)  of  satellites  1  and  3  from  DSST  and  Ji 
model  for  a  800  km,  polar  reference  orbit.  The  satellites  axe  in  a  1  km  circular  formation. 
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Figure  4:  Comparison  of  the  cross-track  separation  (z)  of  satellites  1  and  3  from  DSST  and  Ji  model 
for  a  800  km,  polar  reference  orbit.  The  satellites  are  in  a  1  km  circular  formation. 


Figure  5:  Comparison  of  the  total  separation  ( p )  of  satellites  1  and  3  from  DSST  and  J2  model  for 
a  800  km,  polar  reference  orbit.  The  satellites  are  in  a  1  km  circular  formation. 
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Figure  6:  Comparison  of  the  radial  relative  velocity  (x)  of  satellites  1  and  3  from  DSST  and  J2 
model  for  a  800  km,  polar  reference  orbit.  The  satellites  are  in  a  1  km  circular  formation. 


The  maximum  errors  in  the  J2  model  are  only  a  few  mm/s.  The  errors  shown  here  are  typical  of  all 
the  cases. 

Fig.  10  shows  the  results  for  the  range  rate  between  the  satellites.  The  error  in  the  J2  model 
is  about  30  mm/s  after  four  days.  This  error  is  primarily  caused  by  the  error  in  the  y  component. 
Once  again,  the  error  results  are  typical. 

The  equations  in  the  J2  model  predict  that  there  will  be  no  secular  growth  in  the  relative 
separation  at  critical  inclination  for  satellites  that  have  out-of-plane  separation  caused  only  by 
differences  in  ft.  Figs.  11  through  13  show  that  the  growth  in  relative  separation,  relative  velocity, 
and  range  rate  is  small  according  to  DSST. 


Satellites  1  and  2 


Satellites  1  and  2  primarily  attain  out-of-plane  motion  through  differences  in  inclination.  Con¬ 
sequently,  the  differential  effects  of  J2  are  important.  According  to  the  J2  model  the  largest  change 
in  separation  should  occur  near  an  inclination  of  45°  find  be  mostly  in  the  along-track  component. 
DSST  runs  showed  this  to  be  true.  The  separations  me  compared  in  Figs.  14  and  15.  Fig.  14  shows 
the  large  secular  growth  in  the  along-track  term  which  leads  to  a  large  change  in  the  total  separation. 
Fig.  16  shows  the  range  rate  comparison.  This  is  the  maximum  growth  in  range  rate  for  the  1  km 
circular  formations. 

Fig.  17  shows  the  cross-track  component  for  satellites  1  and  2  for  a  polar  reference  orbit.  The 
growth  in  the  component  is  relatively  small,  but  this  is  the  inclination  at  which  the  maximum  cross¬ 
track  growth  occurs.  Fig.  18  shows  the  relative  velocity.  This  is  the  maximum  relative  velocity  for 
the  1  km  circular  formations. 
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Figure  7:  Comparison  of  the  along-track  relative  velocity  (y)  of  satellites  1  and  3  from  DSST  and 
Ji  model  for  a  800  km  polar  reference  orbit.  The  satellites  are  in  a  1  km  circular  formation. 


Figure  8:  Comparison  of  the  cross-track  relative  velocity  ( z )  of  satellites  1  and  3  from  DSST  and 
model  for  a  800  km,  polar  reference  orbit.  The  satellites  Eire  in  a  1  km  circular  formation. 
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Figure  9:  Comparison  of  the  relative  velocity  («)  of  satellites  1  and  3  from  DSST  and  Ji  model  for 
a  800  km  polar  reference  orbit.  The  satellites  are  in  a  1  km  circular  formation. 


Figure  10:  Comparison  of  the  range  rate  (p)  of  satellites  1  and  3  from  DSST  and  Ji  model  for  a  800 
km,  polar  reference  orbit.  The  satellites  are  in  a  1  km  circular  formation. 
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Figure  13:  Comparison  of  the  range  rate  (p)  of  satellites  1  and  3  from  DSST  and  J2  model  for  a  800 
km,  critically  inclined  reference  orbit.  The  satellites  are  in  a  1  km  circular  formation. 


Time,  hours 

Figure  14:  Comparison  of  the  along-track  separation  ( y )  of  satellites  1  and  2  from  DSST  and  J2 
model  for  a  800  km,  45°  inclination  reference  orbit.  The  satellites  are  in  a  1  km  circular  formation. 
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Figure  15:  Comparison  of  the  total  separation  (p)  of  satellites  1  and  2  from  DSST  and  J2  model  for 
a  800  km,  45°  inclination  reference  orbit.  The  satellites  are  in  a  1  km  circular  formation. 


Time,  hours 

Figure  16:  Comparison  of  the  range  rate  (p)  between  satellites  1  and  2  from  DSST  and  Ji  model 
for  a  800  km,  45°  inclination  reference  orbit.  The  satellites  are  in  a  1  km  circular  formation. 


2180 


Figure  17:  Comparison  of  the  cross-track  separation  ( z )  between  satellites  1  and  2  from  DSST  and 
Ji  model  for  a  800  km,  polar  reference  orbit.  The  satellites  are  in  a  1  km  circular  formation. 


Figure  18:  Comparison  of  the  relative  velocity  («)  between  satellites  1  and  2  from  DSST  and  Ji 
model  for  a  800  km,  polar  reference  orbit.  The  satellites  axe  in  a  1  km  circular  formation. 
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Summary 


The  largest  contributor  to  the  error  in  the  J2  model  compared  to  DSST  was  consistently  the 
along-track  position  component.  Running  a  case  with  J2  only  showed  that  only  a  small  portion 
of  the  error  is  unmodeled  J2  effects.  Most  likely  these  are  higher  order  Jf  effects.  Most  of  the 
along-track  error  is  caused  by  umodeled  solar  radiation  pressure.  This  was  shown  by  a  case  run 
without  solar  radiation  pressure.  The  error  was  much  lower  for  the  case  run  with  SRP  turned  off. 
The  epoch  is  near  the  autumnal  equinox  and  the  orbit  plane  is  nearly  parallel  to  the  sun-earth 
vector.  This  explains  why  the  effect  is  mostly  in  the  along-track  direction.  The  effects  of  SRP  need 
further  investigation  to  assure  that  they  are  real.  Regardless,  the  accuracy  of  the  model  should 
improve  if  ft  is  changed  or  an  epoch  not  so  close  to  equinox  is  used.  A  projected  circular  case  was 
run  at  600  km  altitude.  Unexpectedly,  the  relative  position  errors  were  slightly  smaller  at  the  lower 
altitude.  However,  this  can  be  explained  by  the  fact  that  SRP  was  the  primary  cause  of  the  errors. 
Projected  circular  formations  with  radii  of  100  m  and  10  km  were  also  examined.  All  the  errors  were 
much  larger  for  the  10  km  projected  circular  formation  because  the  greater  separations  translate 
into  greater  differential  perturbations.  However,  the  100  m  projected  circular  formation  still  had 
errors  as  large  as  the  1  km  projected  circular  formations.  In  general,  the  J2  model  was  accurate  over 
four  days  to  within  1  m  in  the  radial  component,  a  few  meters  in  cross-track,  and  less  than  50  m  in 
along-track  and  range.  The  accuracy  in  the  velocity  components  was  better  than  10  mm/s.  Finally, 
the  range  rate  accuracy  was  less  than  1  cm/s.  The  exception  to  these  accuracy  comparisons  was  the 
formation  that  had  a  10  km  nominal  separation.  The  maximum  errors  in  each  component  for  the 
all  the  cases  examined  are  given  in  the  Appendix. 

The  extrema  of  the  relative  motion  from  DSST  is  also  contained  in  a  table  in  the  Appendix. 
As  expected  the  maximum  perturbations  were  for  the  formation  with  the  largest  separation.  The 
formations  at  an  inclination  of  45°  had  the  largest  changes  in  separation  and  range  rate  for  the  1 
km  formations.  In  fact,  the  separation  was  as  low  as  50  m  and  as  large  as  2.56  km  for  the  projected 
circular  formation  at  45°  inclination.  However,  the  polar  orbits  had  the  largest  relative  velocities. 


CONCLUSIONS 


A  model  that  over  severed  days  accurately  describes  the  mean  relative  motion,  relative  velocity, 
and  range  rate  of  satellites  flying  in  formation  has  been  developed.  The  model  includes  the  absolute 
and  differential  effects  of  J2  on  the  dynamics.  The  absolute  effects  of  J2  cause  the  in-plane  and 
out-of-plane  motion  to  have  different  frequencies.  This  causes  a  long  periodic  variation  in  the  sepa¬ 
ration  between  the  satellies.  For  satellites  without  inclination  differences  this  is  the  primary  effect. 
Differential  J2  effects  occur  for  satellites  with  inclination  differences.  These  cause  the  separation 
between  the  satellites  to  grow  secularly.  The  model  has  been  tested  by  comparing  the  results  from 
mean  element  propagations  performed  using  the  Draper  Semianalytic  Satellite  Theory  for  a  number 
of  inclinations,  altitudes,  formation  types,  and  formation  sizes.  The  primary  error  source  not  mod¬ 
eled  was  solar  radiation  pressure  which  caused  an  unmodeled  along-track  drift.  The  model  provides 
a  great  improvement  over  using  the  unperturbed  Hill’s  equations  and  provides  physical  insight  into 
the  errors  in  using  Hill’s  equations. 

In  addition,  the  relative  motion,  relative  velocity,  and  range  rate  expected  between  satellites 
flying  in  formation  has  been  quantified  for  a  number  of  formations.  This  information  will  be  im¬ 
portant  for  the  design  of  relative  motion  sensors  such  as  intersatellite  range  and  range  rate  devices. 
This  will  provide  for  better  relative  motion  and  velocity  estimation  which  is  a  critical  technology  for 
flying  satellites  in  close  formations. 
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There  are  several  areas  of  future  work  related  to  this  effort.  The  bias  in  the  initial  conditions 
'  the  Ji  model  needs  to  be  reflected  in  the  initial  conditions  used  in  DSST.  Also,  the  unusual  SRP 
fects  need  further  investigation  to  see  if  the  effects  are  real.  Jf>  ^3>  and  J4  perturbations  on  ft, 
,  and  M  can  be  added  into  the  equations.  Finally,  the  equations  developed  here  are  valid  for  low 
centricity.  Expanding  the  equations  to  highly  eccentric  orbits  needs  investigation. 
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Table  5:  Maximum  errors  in  J2  model  relative  motion  compared  to  DSST 


Error  in: 

Formation 

X 

[m] 

y 

M 

z 

M 

P 

[m] 

X 

[mm/s] 

y 

[mm/s] 

i 

[mm/s] 

V 

[mm/s] 

P 

[mm/s] 

Circ. 

'  *  =  1° 

Satsl-2 

34.14 

31.69 

1.78 

8.96 

6.96 

Satsl-3 

45.31 

8.18 

43.15 

0.87 

1.90 

8.33 

5.33 

58.55 

i  =  28.5° 

Satsl-2 

0.58 

35.14 

1.60 

33.91 

1.07 

5.50 

2.15 

4.61 

45.87 

Satsl-3 

0.41 

51.12 

3.43 

48.18 

0.55 

1.36 

3.54 

2.82 

72.76 

J2  Only 

Satsl-2 

0.31 

5.48 

1.88 

5.33 

0.50 

4.35 

2.15 

4.23 

5.35 

Satsl-3 

0.62 

12.55 

2.51 

11.22 

0.76 

1.49 

2.90 

2.65 

18.93 

No  SRP 

Satsl-2 

0.61 

8.26 

1.57 

6.74 

1.06 

5.45 

2.15 

4.28 

7.60 

Satsl-3 

0.45 

10.05 

3.43 

9.57 

0.53 

1.08 

3.53 

2.79 

16.12 

i  —  45° 

Satsl-2 

0.50 

33.43 

2.35 

33.12 

0.67 

5.21 

1.89 

4.95 

79.15 

Satsl-3 
i  =  63.435° 

0.34 

41.59 

1.60 

40.44 

0.28 

0.87 

1.78 

1.63 

48.09 

Satsl-2 

0.90 

20.02 

1.71 

18.83 

0.93 

4.99 

3.03 

4.17 

21.61 

Satsl-3 

0.34 

20.11 

1.67 

19.25 

0.22 

0.52 

1.50 

1.15 

20.31 

i  =  90° 

Satsl-2 

0.44 

26.31 

3.67 

25.46 

0.68 

1.00 

5.25 

4.30 

36.56 

Satsl-3 

0.58 

27.10 

2.50 

26.02 

0.54 

1.31 

1.98 

0.84 

30.32 

Proj.  Circ. 

*  =  1° 

Satsl-2 

0.39 

33.97 

9.30 

31.28 

0.72 

1.75 

10.36 

8.32 

61.15 

Satsl-3 

0.63 

47.50 

9.47 

45.10 

0.83 

1.85 

9.64 

6.60 

67.39 

i  =  28.5° 

Satsl-2 

0.62 

33.74 

1.85 

32.55 

1.10 

6.05 

2.48 

5.10 

48.76 

Satsl-3 

0.48 

52.90 

4.00 

49.92 

0.54 

1.28 

4.12 

3.39 

74.37 

Alt=600km 

Satsl-2 

1.54 

30.64 

3.45 

27.53 

2.33 

8.27 

3.96 

5.27 

33.62 

Satsl-3 

1.32 

46.08 

6.40 

41.82 

0.99 

2.63 

6.83 

4.84 

65.55 

Po=10  km 

Satsl-2 

71.19 

135.29 

15.07 

106.51 

33.41 

94.34 

28.41 

60.46 

190.16 

Satsl-3 

17.74 

169.90 

40.21 

159.27 

10.09 

14.07 

39.49 

29.02 

243.81 

Po=100  m 

Satsl-2 

0.44 

40.41 

0.22 

39.07 

0.23 

1.22 

0.25 

0.84 

51.53 

Satsl-3 

0.39 

42.21 

0.42 

40.42 

0.20 

0.77 

0.45 

0.76 

58.71 

i  =  45° 

Satsl-2 

0.61 

33.10 

2.74 

32.82 

0.71 

5.85 

2.21 

5.56 

82.92 

Satsl-3 

0.31 

42.71 

1.84 

41.37 

0.27 

0.81 

2.04 

1.88 

44.95 

i  =  63.435° 

Satsl-2 

1.00 

20.44 

1.98 

19.17 

1.01 

5.69 

3.50 

4.70 

32.14 

Satsl-3 

0.33 

21.01 

1.93 

19.97 

0.25 

0.51 

1.73 

1.39 

19.13 

i  =  90° 

Satsl-2 

0.43 

26.31 

4.26 

25.42 

0.65 

1.05 

6.08 

5.17 

38.63 

Satsl-3 

0.58 

27.26 

2.89 

26.12 

0.54 

1.33 

2.29 

1.07 

28.38 

—  - - - - - - - - - - — — - - - - 
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Table  6:  Extrema  from  DSST  for  different  formations 


Extrema  in: 

Formation: 

Pmin 

[km] 

Pm  a  x 
[km] 

V 

[m/s] 

P 

[m/s] 

Circ. 
i=  1° 

Satsl-2 

0.64 

1.27 

Satsl-3 

0.63 

1.29 

BIS'iis 

t  =  28.5° 

Satsl-2 

0.19 

2.18 

1.22 

1.13 

Satsl-3 

0.72 

1.22 

1.25 

0.52 

i  —  45° 

Satsl-2 

0.11 

2.35 

1.09 

1.08 

Satsl-3 

0.84 

1.15 

1.16 

0.27 

i  =  63.435° 

Satsl-2 

0.14 

2.12 

1.18 

1.13 

Satsl-3 

0.98 

1.02 

1.04 

0.02 

i  =  90° 

Satsl-2 

0.74 

1.25 

1.31 

0.54 

Satsl-3 

0.89 

1.10 

1.12 

0.19 

Proj.  Circ. 
f  =  l° 

Satsl-2 

0.66 

1.36 

1.39 

0.71 

Satsl-3 

0.63 

1.37 

1.40 

0.76 

i  =  28.5° 

Satsl-2 

0.20 

2.38 

1.32 

1.22 

Satsl-3 

0.76 

1.32 

1.33 

0.56 

Alt=600km 

Satsl-2 

0.17 

2.60 

1.38 

1.29 

Satsl-3 

0.72 

1.33 

1.41 

0.65 

pj,=10  km 

Satsl-2 

2.15 

24.19 

13.18 

12.25 

Satsl-3 

7.77 

12.95 

13.29 

5.33 

Po=100  m 

Satsl-2 

0.02 

0.23 

0.13 

0.12 

Satsl-3 

0.05 

0.16 

0.13 

0.09 

O 

ii 

'*» 

Satsl-2 

0.05 

2.56 

2.56 

1.16 

Satsl-3 

0.87 

1.21 

1.25 

0.35 

i  =  63.435° 

Satsl-2 

0.14 

2.29 

1.29 

1.25 

Satsl-3 

0.99 

1.12 

1.16 

0.13 

i  —  90° 

Satsl-2 

0.78 

1.36 

1.41 

0.61 

Satsl-3 

0.93 

1.17 

1.21 

0.25 
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